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METHODS OF EXTRACTING SPICE PARAMETERS, PERFORMING A 
SPICE CALCULATION, AND PERFORMING DEVICE SIMULATION FOR 
A PARTIALLY-DEPLETED SOI MOSFET 

BACKGROUND OF THE INVENTION 

[0001] Field of The Invention 

[0002] The present invention relates to a MOSFET formed on an 
SOI (Silicon On Insulator) substrate, and in particular to a method of 
detecting SPICE (Simulation Program With Integrated Circuit Emphasis) 
parameters, performing a SPICE calculation, and performing device 
simulation for a partially depleted SOI MOSFET where the body is only 
partially depleted. 

[0003] Related Art 

[0004] An SOI MOSFET is a MOSFET constructed in a single 
silicon crystal formed on an insulator film called an SOI (Silicon On Insulator) 
substrate. In particular, in a partially depleted SOI MOSFET with a floating 
body effect, each component of the currents that flow into a body current and 
causes changes in the amount of body charge (such currents are collectively 
referred to as the "body current") is extremely small, which means that the 
fluctuations in the body charge are also extremely small. This means that 
during circuit operation, capacitive couphng occurs between the body and 
each of the nodes, i.e., the gate, the source, the drain, and the back gate. It is 
also known that with methods where a MOSFET is operated by providing a 
terminal on the body and fixing the potential, capacitive coupling occurs in the 
same way when the body or body terminal resistance and circuit operating 
fi:equency are both sufficiently high. 

[0005] In this way, in an SOI MOSFET where capacitive coupling 
occurs in the body, during actual circuit operation, the total body charge 
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affects not only the body potential but also the effective drain current 
characteristics. Naturally, when the effective draiin current characteristics 
are estimated by performing device simulation and a SPICE calciilation, it is 
necessary to calculate the total body charge during circuit operation. 

[0006] An article by A. Wei et al. in IEEE Trans. Electron Devices, 
voL45, pp. 430-438, 1998 discloses a method of calcxilating, through device 
simulation, the total body charge of an SOI MOSFET during circuit operation. 

[0007] In SPICE calculations of an SOI MOSFET, it is normally 
thought necessary to prepare models that reproduce the body current, that is, 
a diode current, an impact ionization current, and a tunnel current resulting 
in gate induced drain leakage (GIDL), and the capacitance-voltage 
characteristics (or charge -voltage characteristics) of the body. 

[0008] For example, kinks that appear in the SOI MOSFET static 
characteristics (the output characteristics and the sub-threshold 
characteristics) are understood as a phenomenon where the potential of the 
body increases due to the impact ionization current- In order to reproduce 
this kink phenomenon in simulations, models of the current-voltage 
characteristics for the body current and body capacitance -voltage 
characteristics have been implemented in the standard models that are 
presently available (the SOI version of BSIM3 caUed "BSIMPD" by UC 
Berkeley and UFSOI by U. Florida, etc.). 

[0009] However, with regard to models of the current-voltage and 
capacitance -voltage characteristics for the body (substrate) that are not so 
important in conventional bulk-MOSFET techniques, the need to perform a 
procedure for estimating and calculating coefQcients for specifying the details 
of these models in keeping with actual device characteristics presents a 
considerable burden. 

[OOlO] Since the body current that is an origin of the fluctuation in 
the total body charge during actual circuit operation is extremely small, after 
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sufficient repetition of switching operations, the total charge of the body 
converges on a state where the charge fluctuates slightly about a certain value. 
During circuit simiilation, it is necessary to calculate the convergence value 
for the total body charge in advance and to set this calculated value as the 
initial value for simulation. In addition, since the absolute value of the body 
current is an order of magnitude smaller than that of the driving current, the 
calcvdation of the total body charge and moreover the potential need to be 
calculated with high accuracy such that there are inevitable compromises 
between the time taken for the calculation (calculation time) and the degree of 
reliability of the convergence. 

[001 1] When performing circuit analysis on a partially depleted SOI 
MOSFET, it is necessary to estimate and calculate coefficients that determine 
the details of models of the current-voltage and capacitance-voltage 
characteristics for the body in order to meet actual characteristics, and to 
perform a procedure for calculating in advance the convergence value of the 
total body charge during circuit operation. That is, the procedure for 
extracting the SPICE parameters and performing a SPICE calculation is 
complex. While the calcvdations of the total charge and moreover the 
potential of the body need to be performed with high accuracy, the calculation 
time and convergence unavoidably have to be prioritized, so a compromise is 
made to the calculation acciuracy. 

[0012] The present invention was conceived in view of the situation 
described above, and attempts to increase the efficiency of the procedures for 
SPICE parameter extraction, SPICE calculations, and device simulation for a 
partially depleted SOI MOSFET. 

SUMMARY 

[0013] A method of extracting SPICE parameters and performing a 
SPICE calculation according to the present invention relates to a partially 
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depleted SOI MOSFET formed on an SOI substrate, wherein a body cxirrent is 
not modeled and is estimated at zero, and information regarding a steady 
state dviring circuit operation is treated as macro parsuneters that encompass, 
on a macro level, amounts that are originally to be found, thereby omitting an 
extraction of parameters relating to the body current and a calculation that 
finds the steady state during circuit operation. 

[0014] A method of performing a SPICE calculation for a partially 
depleted SOI MOSFET according to the present invention relates to a 
partially depleted SOI MOSFET formed on an SOI substrate, wherein a body 
current is estimated at zero and information regarding a steady state during 
circuit operation is treated as macro parameters that encompass, on a macro 
level, amounts that are originally to be found, thereby omitting a procedure 
for a calculation that finds the steady state during circuit operation. 

[0015] According to the respective methods of extracting SPICE 
parameters and performing a SPICE calculation for a partially depleted SOI 
MOSFET in accordance with the present invention described above, 
information in the steady state during circuit operation is intentionally 
treated as macro parameters instead of being found from the body current 
characteristics according to customary principles. By doing so, it is possible 
to omit an extraction of the SPICE parameters relating to the body current 
and a calculation that finds the steady state during circuit operation. This 
can be regarded as a reasonable means in view of the tradeoff between the 
calculation accuracy and the calculation time and convergence, which 
originates from the absolute value of the body current being small. 

[0016] A method of performing a SPICE calculation and device 
simulation according to the present invention relates to a partially depleted 
SOI MOSFET formed on an SOI substrate, wherein an input signal cycle and 
an output load capacitance or an output load resistance are all increased by a 
same ratio "X", an initial value '^Y" of a total body charge at which there is zero 



-4- 



9319S-000577 



fluctuation in the total body charge after one input signal cycle is found in 
advance, and a value of Y for when X=l is extrapolated by expressing Y as a 
function of X, thereby finding a steady state dxuing circuit operation of a 
partially depleted SOI MOSFET. 

[0017] A method of performing a SPICE calculation and device 
simulation according to the present invention relates to a partially depleted 
SOI MOSFET formed on an SOI substrate, wherein a cycle of trial pulses that 
imitate actual circuit operation is increased by "X" times to extend the trial 
pulses with similarity waveforms, an initial value 'T' of a total body charge at 
which there is zero fluctuation in the total body charge after one trial pulse 
cycle is found in advance, and a value of Y for when X=l is extrapolated by 
expressing Y as a function of X, thereby finding a steady state during circuit 
operation of a partially depleted SOI MOSFET. 

[0018] According to the respective methods of performing a SPICE 
calcxilation and device simulation for a partially depleted SOI MOSFET in 
accordance with the present invention described above, only the cycle is 
increased by X times without changing the history pattern of the oscillations 
in the bias conditions. It is known that the initial value 'Y" of the total body 
charge at which there is zero fluctuation in the total body charge after one 
cycle of the input signal is highly dependent on only the history pattern of the 
oscillations in the bias conditions. That is, the X dependency of the value Y is 
extremely small, which suggests that it will be extremely easy to extrapolate 
the value of Y for when X=l. Also, increasing the cycle of the oscillations in 
the bias conditions by X times means that the accuracy of the calculation of 
the total body charge is also increased by X times. Accordingly, the accuracy 
of the calculation can be increased without increasing the calculation time and 
the load of the calculation procedure. 

[0019] A method of performing a SPICE calculation according to the 
present invention relates to a partially depleted SOI MOSFET formed on an 
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SOI substrate, wherein both special values set with lower than a normal 
magnitude for calcxxlation for a body, and values set with a normal magnitude 
that is used in a calculation for nodes aside from the body are utilized as error 
tolerance parameters used as a convergence criterion. 

[0020] According to the method of performing a SPICE calculation 
for a partially depleted SOI MOSFET in accordance with the present 
invention described above, it is possible to perform a high accuracy calculation 
only for the body. By doing so, it is possible to cope with a situation when the 
body current and moreover the fluctuations in total body charge are extremely 
small and when it is necessary to perform the calculation of the total body 
charge and moreover the body potential with high accuracy. At the same 
time, since the calculation of the charging and discharging by the driving 
current at nodes aside from the body can be performed with normal accuracy, 
this method is more efficient than calculating all of the nodes, including the 
body, with high accuracy. 

[0021] A method of performing a SPICE calcvdation for a partially 
depleted SOI MOSFET according to the present invention relates to a 
partially depleted SOI MOSFET formed on an SOI substrate, wherein for a 
partially depleted SOI MOSFET which is provided with a body terminal and 
whose effective channel width is known, when body potential dependence of 
drain current characteristics in the steady state is measured and modeled, a 
body potential — a value at which a drain current per unit effective channel 
width matches that in a steady state of a floating body partially depleted SOI 
MOSFET whose effective channel width is known — is found in advance 
directly or by interpolation/extrapolation, and the body potential is 
tabiilarized and used in analysis of the steady state for the floating body 
partiaUy depleted SOI MOSFET. 

[0022] According to the method of performing a SPICE calculation 
for a partially depleted SOI MOSFET in accordance with the present 
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invention described above, it is possible to find the body potential in the 
steady state even in cases where the body current is not modeled. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0023] FIG. 1 is a flowchart showing a method of performing a 
SPICE calculation for a partially depleted SOI MOSFET according to the first 
embodiment of the present invention. 

[0024] FIG. 2 is a flowchart showing a method of performing a 
SPICE calculation and device simvilation for a partiaUy depleted SOI 
MOSFET according to the second embodiment of the present invention. 

[0026] FIG. 3 is a flowchart showing a method of performing a 
SPICE calculation for a partially depleted SOI MOSFET according to the 
third embodiment of the present invention. 

[0026] FIG. 4 is a cross-sectional drawing showing a partially 
depleted SOI MOSFET. 

DETAILED DESCRIPTION 

[0027] FIG. 1 is a flowchart showing a method of extracting SPICE 
parameters and performing a SPICE calculation for a partially depleted SOI 
MOSFET according to a first embodiment of the present invention. 

[0028] FIG. 4 is a simplified cross-sectional drawing showing a 
partially depleted SOI MOSFET. A single siUcon crystal is formed on an Si02 
film 41 called a "buried oxide film" or "box layer", and it comprises a p-type 
body 42 and an n-type source/drain region 44 and a channel 43 for example. 
A gate electrode 46 is constructed on the channel 43 via a gate oxide film 45. 
Side walls (spacers) 47 are formed on both sides of the gate electrode 46 after 
forming an N-type extension region 441 with a lower concentration than the 
source/drain region 44. The above describes a floating body partially 
depleted SOI MOSFET. The broken line in FIG. 4 shows the construction of 
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a partially depleted SOI MOSFET in the case where a body terminal 51, i.e., a 
terminal that extends out of the body 42, is provided. 

[0029] In the present invention, a model is produced for the 
electrical characteristics of a partially depleted SOI MOSFET in a steady 
state during circuit operation. The expression "steady state dxiring circuit 
operation" refers to a state where the body charge does not change even after 
one cycle of the circuit operation. If Qb is the body charge, the following 
expression can be defined- 

Qb(t+tp) - Qb(t) = 0 ... (1). 

[0030] Here, Qb(t) represents the body charge at an arbitrary time t. 
Also, tp represents one cycle of the circuit operation. 

[0031] In FIG. 1, first predetermined parameters, such as extracted 
parameters that describe the electrical characteristics of a partially depleted 
SOI MOSFET, a condition regarding the supply voltages and input signals 
that are based on a predetermined circuit, and error tolerance parameters 
that are used when judging convergence for the SPICE calculation, are set in 
conformity with a corresponding descriptive style (process Sll). However, at 
this time, the fluctuations in the body charge that are observed in one cycle of 
actual circuit operation are thought to be sufficiently small, so that 
parameters describing the body current characteristics are not extracted and 
the body current is estimated at zero (process Slla). The parameters that 
relate to the body current include, as examples, parameters relating to a diode 
current, parameters relating to a parasitic BJT (Bipolar Jimction Transistor) 
current, parameters relating to an impact-ionization current, and parameters 
relating to a GIDL current. 

[0032] Then, in place of the parameters that describe the body 
current characteristics, information in the steady state during circuit 
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operation that is normally found from the body current characteristics, which 
is to say, the body charge and oscillation in the body potential, is treated as 
macro parameters that encompass information on the body current 
characteristics (process SI lb). After setting parameters including such 
macro parameters, a SPICE calculation for transient analysis is performed 
(process S12). It is noted that each macro parameter should be determined 
from corresponding input information for the SPICE calcvilation by using 
interpolation and/or extrapolation. 

[0033] With the method of the present embodiment described above, 
information regarding the steady state dxiring circuit operation is 
intentionally handled as macro parameters instead of finding such 
information from the body current characteristics based on elementary 
principles- By doing so, it is possible to omit the extraction of the SPICE 
parameters relating to the body current and a calculation that finds the 
steady state during circuit operation. This can be regarded as a reasonable 
means in view of the tradeoff between the calculation accuracy and the 
calculation time and convergence, which originates from the absolute value of 
the body current being small. 

[0034] As one example, when calculating the propagation delay of a 
CMOS inverter circxiit that is composed of a partially depleted SOI MOSFET, 
the fluctuation in the body charge observed in one cycle of actual circuit 
operation is thought to be sufficiently small so that transient analysis can be 
performed with macro parameters that estimate the body current at zero. 
When performing transient analysis for a normal short period, if the error 
tolerance parameters are set high, it is possible to easily prioritize both the 
calculation time and convergence. 

[0036] FIG. 2 is a flowchart showing a method of performing a 
SPICE calcidation and device simulation of a partially depleted SOI MOSFET 
according to a second embodiment of the present invention. Basically, in 
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keeping with the article by A. Wei et al. in IEEE Trans. Electron Devices, 
vol.45, pp. 430-438, 1998, a plurality of internal device states with different 
total amounts of body charge are prepared as the initial states for analysis 
(process S21). After this, input signals or trial pulses that imitate the actual 
circuit operation are apphed to corresponding nodes and the initial body 
charge whose change after one cycle is zero (which is determined as the 
effective body charge though the following procedure) is used in the method of 
the present embodiment. 

[0036] That is, as shown in process S22, for the input signals or trial 
pulses that imitate the actual circuit operation, the cycle is increased by "X" 
times to produce elongated waveforms with similarity forms. These 
waveforms are apphed to the corresponding nodes. By doing so, the initial 
value for the total body charge at which the fluctuation after one cycle 
becomes zero is foimd (process S23). After this, Y is expressed as a function 
of X, and the value of Y when X=l is extrapolated (process S24). In process 
S24, the extrapolated value of Y is determined to be the effective value of the 
total body charge during actual circuit operation. 

[0037] According to the method of the present embodiment described 
above, during a SPICE calculation and device simulation for a partially 
depleted SOI MOSFET, the cycle of the input signals during actual circuit 
operation or of trial pulses that imitate such signals is increased by X times to 
form an elongated waveform for which an initial value 'IT' of the total body 
charge is conserved. After this, the analysis result for the case where X=l, 
i.e., when the waveform is not elongated, is found by extrapolation. When 
the waveform is elongated, only the cycle is increased by X times, and the 
history pattern of the oscillations in the bias conditions is not changed. It is 
known that the initial value Y of the total body charge at which there is zero 
fluctuation in the total body charge after one cycle of the input signal is highly 
dependent on only the history pattern of the oscillations in the bias conditions. 
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That is, the X dependency of the value Y is extremely small, which suggests 
that it will be extremely easy to extrapolate the value of Y for a case where 
X=l. Also, increasing the cycle of the oscillations in the bias conditions by X 
times means that the accuracy of the calculation of the total body charge is 
also increased by X times. Accordingly, the accxiracy of the calculation can be 
increased without increasing the calculation time and the load of the 
calcxilation procedure. 

[0038] Also, the accuracy of the calculation may be improved using 
the following method in place of the method described above. Both special 
values set with lower than a normal magnitude for calculation for the body, 
and values set with a normal magnitude that is used in a calcidation for the 
nodes aside from the body are utiHzed as error tolerance parameters used as a 
convergence criterion of the SPICE calculation. According to this method, it 
is possible to perform the high accuracy calculation only for the body. By 
doing so, it is possible to cope with a situation when the body current and 
moreover the fluctuations in total body charge are extremely small and when 
it is necessary to perform the calculation of the total body charge and 
moreover the body potential with high accuracy. At the same time, since the 
calculation of the charging and discharging by the driving ciurent at nodes 
aside from the body can be performed with normal acciu-acy, it can be said 
that this method is more efficient than performing a calculation based on all of 
the nodes, including the body, and the calculation is performed with high 
accuracy. 

[0039] FIG. 3 is a flowchart showing the method of performing a 
SPICE calculation for a partially depleted SOI MOSFET according to a third 
embodiment of the present invention. In the first embodiment described 
above, in place of parameters that describe the body current characteristics, 
information regarding the steady state during circuit operation that is 
normgdly found from the body current characteristics, that is, the body charge 
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and body potential oscillation, is treated as macro parameters that encompass 
information on the body current characteristics. On the other hand, in this 
third embodiment, a method is described where in place of parameters that 
describe the body current characteristics, the body potential not during the 
circuit operation but dining what is normally referred to as the "steady state" 
is tabularized and used in a SPICE calcxdation. 

[0040] The drain current characteristics are measured for a T-gate, 
H-gate or similar type of partially depleted SOI MOSFET for which the 
effective channel width is known and which is provided with a body terminal. 
Alternatively, if a model has been prepared for the same type of SOI MOSFET, 
such model is used (process S3l). Then, the body potential - a value at which 
a drain current per unit effective channel width matches that in the steady 
state of a floating body partially depleted SOI MOSFET whose effective 
channel width is known — is found directly or through 
interpolation/extrapolation (process S32). In addition, the calculated body 
potential is tabvilarized and used in the analysis of the steady state for a 
floating body partially depleted SOI MOSFET (process S33). 

[0041] With the method of the above embodiment, it is possible to 
calcxdate the body potential in the steady state, even when the body current 
has not been modeled as shown in the first embodiment. Accordingly, an 
increase in calculation efficiency can be expected. 

[0042] According to the method of the present embodiment, in 
transient analysis that takes the body current into account, the electrical 
characteristics diiring circuit operation are heavily dependent on the error 
tolerance parameters used during the convergence determination of a SPICE 
calculation. Also, a measure that sets the error tolerance parameters low to 
increase the accuracy of a SPICE calciilation sacrifices calculation time and 
convergence, and is not necessarily effective. In contrast, by estimating the 
body current at zero and performing transient analysis with macro 
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parameters, the error tolerance parameters are set high and the calculation 
time and convergence are prioritized. The total body charge is determined so 
as to be appropriate for the propagation delay resulting from actual 
measurement, and a suitable oscillation pattern for the body potential is 
determined. By doing so, SPICE calculations can be made more e£Eicient. It 
shovdd be noted that for a partially depleted SOI MOSFET provided with a 
terminal in the body, in cases where the body resistance and circuit operating 
frequency are both sufficiently high for transient capacitive coupling to occur, 
it is possible to increase efficiency with exactly the same procedure. 

[0043] Also, by increasing the cycle of the input signals or trial 
pxdses that imitate the actual circuit operation, invariance in the total charge 
of the body after one cycle of circuit operation is found and invariance for a 
case when the cycle is thereafter returned to normal is found by extrapolation. 
By doing so, accuracy is improved. Alternatively, both special values set 
with lower than a normal magnitude for calculation for the body, and values 
set with a normal magnitude that are used in the calcidation for the nodes 
aside from the body may be utiUzed as error tolerance parameters used as a 
convergence criterion of a SPICE calculation. It is also possible to calculate 
the body potential in the steady state in advance and to tabularize and use the 
resvdts in a SPICE calculation for the steady state. 

[0044] As described above, according to the present invention, 
instead of finding information regarding the steady state during circuit 
operation from the body current characteristics based on fundamental 
principles, the body current is estimated at zero, and the same information is 
treated as macro parameters. By doing so, it is possible to omit the 
extraction of SPICE parameters relating to the body cxurent characteristics 
and a calcidation that finds the steady state during circxiit operation. In 
view of the tradeoff between the calcvdation accuracy and the calculation time 
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and convergence, which originates from the absolute value of the body current 
being small, this means can be regarded as reasonable. 

[0045] Also, by increasing the cycle of the input signals or trial 
pulses that imitate the actual circuit operation, invariance in the total charge 
of the body after one cycle of circuit operation is found and invariance for a 
case when the cycle is thereafter returned to normal is found by extrapolation. 
By doing so, accuracy is improved. Alternatively, both special values set for 
the calcxilation for the body that is with lower than a normal magnitude and 
values set with a normal magnitude that are used in the C£dculation for the 
nodes aside from the body may be utiUzed as the error tolerance parameters 
used as a convergence criterion of a SPICE calcxilation. It is also possible to 
calculate the body potential in the steady state in advance and to tabxilarize 
and use the results in a SPICE calculation of the steady state. As a result, it 
is possible to improve the efficiency of the procedures for SPICE parameter 
extraction, SPICE calculation, and device simulation for a partially depleted 
SOI MOSFET. 

[0046] The entire disclosure of Japanese Patent Application No. 
2002-345823 filed November 28, 2002 is incorporated by reference. 
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